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ABSTRACT 

We present a new derivation of the CORS Baade- Wesselink method in the Walraven 
photometric system. We solved the complete Baade- Wesselink equation by calibrating 
the surface brightness function with a recent grid of atmosphere models. The new 
approach was adopted to estimate the mean radii of a sample of Galactic Cepheids for 
which are available precise light curves in the Walraven bands. Current radii agree, 
within the errors, quite well with Cepheid radii based on recent optical and near- 
infrared interferometric measurements. 

We also tested the impact of the projection factor on the Period-Radius relation 
using two different values {p — 1.36, p = 1.27) that bracket the estimates available in 
the literature. We found that the agreement of our Period-Radius relation with similar 
empirical and theoretical Period-Radius relations in the recent literature, improves by 
changing the projection factor from p = 1.36 to p = 1.27. Our Period-Radius relation 
is logi? = (0.75 ±0.03) log P+ (1.10 ±0.03), with a rms = 0.03 dex. 

Thanks to accurate estimates of the effective temperature of the selected Cepheids, 
we also derived the Period-Luminosity relation in the V band and we found My = 
(-2.78 ± 0.11) log P + (-1.42 ± 0.11) with rms = 0.13 mag, for p = 1.27. It agrees 
quite well with recent results in the literature, while the relation for p — 1.36 deviates 
by more than 2a. We conclude that, even taking into account the intrinsic dispersion 
of the obtained Period-Luminosity relations, that is roughly of the same order of 
magnitude as the effect of the projection factor, the results of this paper seem to 
favour the value p = 1.27. 

Key words: Cepheids - Baade- Wesselink radii - Period-Radius relation - Walraven 
photometry. 



1 INTRODUCTION 

Classical Cepheids are fundamental standard candles 
to determine both Galactic and extragalactic distances. 
They obey well defined optical and Near-Infrared 
Period-Luminosity and Period-Luminosity-Colour rela- 
tions. Thanks to the unprecedented spatial resolution of the 
Hubble Space Telescope, Cepheids can be used as primary 
distance indicators up to the Virgo cluster (~ 25 Mpc). This 
means that they can also be adopted to provide robust ab- 
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solute calibrations of secondary distance indicators, and in 
turn to estimate the value of the Hubble constant (see e.g. 
iFreedman et aLlbOOlh . 

The accuracy reached in the determination of cosmo- 
logical parameters using the Cepheid distance scale re- 
lies on the knowledge of the physical mechanisms that 
govern their radial oscillations and on their dependence 
on intrinsic paramete rs (luminosity, mass, radius, chemica l 
comp o sition; see e.g. Bono. Marconi fc Stellingwerj (Il999l, 



l2OO0h : iBono. Castellani fc Marconil (|2002l ): iKeller fc Wood 
1 20061 ) ■ and references therein). In particular, the estimate 
of Cepheid radii and of the effective temperature (using 
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mean intrinsic colours) allows us to constrain their intrin- 
sic luminosity. Furthermore, the Cepheid mass can be esti- 
mated by adopting a Period-Mass-Radius relation (see e.g. 
iBono et a l. 2001b, and references therein), and in turn to 
provide independent constraints on the mass-luminosity re- 
lations. 

The methods currently used to derive Cepheid radii 
from photometric and spectroscopic data (radial veloc- 
ity) are ba s ed on the classical Baade-Wesselink technique 
l|Wesselinkl 19461'). Among them , the surface brightness 
technique (iBarnes fc Evand 1 19761 ) and the CORS method 
l|Caccin et al.lll98lh rely on solid physical bases. The reader 
interested in a detailed discussion concer ning the differen t 
Baade-Wesselink methods is referred to iGautschvl l|l987t ). 
The use of interferometric methods and the measurement 
of parallaxes provide a model-independent method for ra- 
dius measurements. However, this approach has only been 
applied to a very limited sample of bright nearby Cepheids 
jNordgren et al. 200Gl : lLane et al ] [200 2l: lKervella et al ]|2004 
iDavis et al.l l2009l ). On the the other hand, the Baade- 
Wesselink methods can be applied to large samples of 
Cepheids and it has been rec ently extended to Magellanic 
Cepheids ()Storm et all l2005l l. However, the precision of 
Baade-Wesselink distances strongly depends on the quality 
of photometric and radial velocity data. 

In a previous paper (|Ruoppo et al.ll2004l ). we adapted 
the CORS method to derive Cepheid radii observed in the 
Stromgren system. However the model atmospheres adopted 
in the quoted investigation to estimate the surface brightness 
did not allow us to include Cepheids with periods longer 
than 12-13 days. The cutoff in the long period range limited 
the accuracy of the Period-Radius relation. 

Recently, IPediceUi et al] (|2009l ) provided new calibra- 
tions b a sed o n a new release of model atmospheres by 
ICastellil (| 19991 ) for the Walrave n p h otom etry of Galactic 
Cepheids originally collected by IPell l| 1978^ 1 . The Walraven 
system includes five bands, (V, B, L, U, W), designed to 
measure the features of hydrogen spectrum. Three of them 
have the centra l wavelength in the ultraviolet region (see 
iLub fc Pelll 19771 ; [Pel & Lub 2007, for more details about this 
photometric system) . In particular the central wavelength of 
the Walraven filters are the following: \v = 5405A, As = 
4280A, Xl = 3825A, \u = 3630A, Aw = 3240A. 

In this investigation we show that this new calibration 
of the VBLUW Cepheid time series data bv lPedicelli et ahl 
l|2009[ ) allows us to derive precise stellar parameters (surface 
gravity and effective temperature) at fixed chemical compo- 
sition and to calibrate with high accuracy the surface bright- 
ness function, a fundamental ingredient for the application 
of the Baade-Wesselink method. This opportunity stimu- 
lated us to apply the CORS ve rsion of the Baade-Wes selink 
technique as formulated by lOnnembo et al.l (|l985l ) and 
IRuoppo et al.l (|2004h . to derive reliable radii for a sample 
of Cepheids having periods from a few days up to ~ 40 days 
and to obtain more accurate Period-Radius and Period- 
Luminosity relations. 

The layout of the paper is the following: the original 
CORS method is introduced in Sec. (2] The grids of atmo- 
sphere models and the calibration of the surface brightness 
in the Walraven system are described in Sec.O In Sec.|4]we 
describe the adopted photometric and spectroscopic data, 
while the estimate of the Cepheids radii and of the Period- 



Radius relation are discussed in Sees. [5] and [G] respectively. 
In Sec. [7] we discuss the derivation of the Period-Luminosity 
relation. Finally, Sec. [S] summarizes the main results of this 
investigation. 



2 THE CORS METHOD 



The original CORS method jCaccin et al.|[l981 ) is a variant 
of the classical Baade-Wesselink technique ( Wesselinl3ll946l ) 
and was developed to derive the radii of pulsating stars. It 
relies on the surface brightness function: 



Sv = rnv + 5 log a 



(1) 



where my is the apparent visual magnitude and a is the 
angular diameter of the star. For a variable star Eq. ([T} can 
be written for each phase (0) along the pulsational cycle. By 
differentiating the above equation with respect to the phase, 
by multiplying the result with a colour index (Cij, where i,j 
are two generic photometric bands) and by integrating over 
the entire cycle, one can write: 



ln|i?o-pP / v{(i)')d(i}'^C[jd<t)- 



B + AB = 



(2) 



where q = j^^q i P is the period, v is the radial velocity 
and p is the radial velocity projection factor. The other two 
terms, B and AB, are the following: 



B 



Cij{(j>)m'y{<j>)d(j) 



AB= I aj{(t>)Sv{(t))d(j, . 

10 



(3) 
(4) 



The projection factor p correlates radial and pulsation 
velocity, i.e. R'{4>) = —p ■ P ■ v{(f>)- The uncertainty affecting 
the estimates of this parameter is still the main source of 
systematic errors in the Baade-Wesselink method. We shall 
discuss the p value in Sec. 15.11 

By solving Eq. ([2]) we estimate the value of the radius 
-Ro at an arbitrary phase (po, while the radius at any phase (j> 
along the pulsational cycle can be evaluated by integrating 
the radial velocity curve between (po and (p. The term B 
present in Eq. © can be easily determined using colour and 
magnitude (B), while the term AB is not directly connected 
to observational data, but approximates the area of the loop 
performed by variable stars in the colour-colour plane. 



Typically the AB term has a small value (10 



10" 



lOnnembo et all 1 19851 ) and in the original Ba ade-Wesselink 
method it is neglected fsee lCaccin et S][i983). However, the 
Cepheid radii estimated by including the AB term in the 
Eq. ([2]) are more accurate than those b ased on the original 
Baade-Wesselink method (jSoUazzo et al. 1981; Ripcpi et al] 
1 1991 12OOOI ; IRuoppo et al.ll2004l ). 

In the next section we describe how to calibrate the 
surface brightness function to estimate the term AB from 
Eq. ©. 
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3 THE MODIFIED CORS METHOD 

3.1 The grids of atmosphere models 

In order to calibrate the surface brightness function, we as- 
sume the validity of the quasi-static approximation (QSAn 
for the Cepheid atmosphere (see also lOnnembo et al.ll 19851) . 
Under this hypothesis, any photometric quantity can be ex- 
pressed as a function of effective temperature and gravity 
{T^ff, Qeff)- Then, we can write: 

Sv ^ Sv(Teff, geff) (5) 

for the surface brightness, and: 

Cij = Cij(Teff, Qeff) (6) 
Ckt = C'kl{Teff, Qeff) 

for two intrinsic colour indices. If the determinant of the 
Jacobian J [Cij, Cki\Teff, Qeff) the last two equations 
can be inverted and we obtain: 

Teff = Teffidj, Ckl) (7) 
9eff = QeffiCij, Ckl) 

and we can also write the surface brightness as a function 
of the two colour indices: 

Sv ~ Sv(Cij, Ckl) ■ (8) 

In general, the invertibility condition is not valid 
over the entire parameter space, since the same pair of 
colours trace different pairs of gravity and temperature 
(geff, Teff). Figure [T] shows the colour deg eneracy present 
in the grid of atmosphere model provided by CastcUi ( 1999|) 
( http:/ /wwwuser. oat. ts. astro, it /castelli/ color s/vbluw. html ) , 
for the Walraven two-colour diagram {V - B) vs {U - 

However, a local solution can be found, provided that 
an appropriate choice of the colours dj and Cki and of the 
range in colour is made. By using the quoted colours {V — B) 
and {U — W), we succeeded in inverting Eq. ((6]) for the range 
of parameters typical of Cepheids, i.e. 0.0 ^ logge// ^ 3.0 
dex and 4500 < Te// 7250 K. 

To verify that real Cepheids are located inside the the- 
oretical grid, we plotted the loops, obtained by fitting the 
observed colours of selected Cepheids, in the quoted colour- 
colour plane. The result for three Cepheids with short-, 
medium- and long-period, is shown in Fig. (2] A glance at 
the data plotted in this figure shows that the three loops 
performed by actual Cepheids are located inside the limits 
of t he predicted grid. 

iLub fc Pell l|l977l ) and |Pei (|l978l ) pointed out that the 
flux in the W band was measured with sufficient accuracy 
only for the brightest Cepheids due to the intrinsic faintness 

^ Tlie quasi-static approximation assumes that the atmosphere 
of the pulsating star can be described, at any time, by a clas- 
sical hydrostatic, plane parallel model in radiative/convective 
equilibrium and in local thermodynamic equilibrium, identified 
by the effective temperature (Te//) a-ud by the effective gravity 

{9eff - -rS- + dir)- 

We remember that in the Walraven photometric system the 
magnitude in a generic band X is defined without the —2.5 fac- 
tor: mx = log -Fx- As a consequence, the colour is defined as 
the difference between the magnitude at longer and at shorter 
wavelength. 
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Figure 1. The grids of model atmosphere plotted in the plane 
of iy — B), [U — W) colours: lines of constant temperature and 
constant gravity axe plotted. 



of Cepheid fluxes in the near-ultraviolet, but also because 
the telescope mirror reflectivity in this band was deterio- 
rated at the epoch of observations. Therefore, we decided to 
perform a further test by using the pair of colours (V — B) 
and {L — U), even if the curves of constant temperature 
and gravity are not orthogonal as in the plane with the 
{U — W) colour. Accurate solutions can be obtained in the 
range 5000 ^ T^ff < 7500 K and 0.0 < logpe// < 4.00 dex. 
However, the loops of Cepheids with periods longer than 
11 days fall outside the grid of models. This means that the 
structural parameters of Cepheids and the surface brightness 
function are less accurate using this colour-colour plane. 
Nevertheless, the ratios of the Cepheid radii with P < 11 
days, estimated using the two different colour-colour planes, 
show that the differences are smaller than 10%. On the basis 
of this evidence, we decided to focus our attention on the 
{V — B) — {U — W) plane (see Appendix [B]for more details) . 



3.2 Derivation of Sv 

The availability of a region across the theoretical grid cov- 
ering the locus occupied by the Cepheid colour-colour loops 
and the validity of the invertibility condition shown in the 
above section provide us with an instrument to write ef- 
fective temperature and gravity as a function of the two 
colour indices and, flnally to calibrate the surface bright- 
ness through Eq. ((8]). 

As a flrst step, we have fitted the theoretical grid shown 
in Fig. [2] by means of polynomials, in order to express tem- 
perature and gravity as a function of {V — B) and {U — W) 
colours, according to Eq. ((7| (see Appendix |X] for de- 
tails about the fitted equations). The surfaces logTe// = 
logT,ff{V-B, U-W) undlogg.ff = log g,ff{V~B, U~ 
W) are shown in Fig. [S] |4l The rms around the fitted equa- 
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V-B 

Figure 2. Colour-colour loops for three selected Cepheids plotted 
onto the theoretical grids. The three Cepheids are SS Set (P = 3.7 
days), RZ Vel (P = 20.4 days) and 1 Car (P = 35.5 days). 




Figure 3. The surface representing logT^^^ as a function of 
iy — B) and {U — W), obtained from the polyn omial fit of the 
atmosphere models (crosses) from lCastellil lll999l) . 



tions amount to 0.0033 dex and 0.15 dex for XogT^ff and 
logSe//, respectively. 

Following the Walraven convention for magnitudes, the 
surface brightness can be derived in the following way: 

Sv = const. +A\ogT,ff - BC{T,ff, g,ff) (9) 

where BC is the Bolometric Correction calculated as 
a function of T^ff and g^ff through a polynomial fit (see 
Appendix A plot of the fit of BC as a function of 

log Tef f and log g^f / is shown in Fig. [5] 

The procedure described here allows us to calculate the 
term AB in the complete CORS formulation (Eq. [2]). 




u-w 



Figure 4. The same as Fig.[3l but for log g^ff. 




Figure 5. The surface representing the Bolometric Correc- 
tion, BC, as a function of logT^^^ and log g^j^j^, obtain ed from 
the polynomial fit of atmosphere models (crosses) from ICastellil 
(.1999.) . 

3.3 Test with spectroscopic and interferometric 
data. 

In the previous section we have described the procedure fol- 
lowed to derive the physical parameters of Cepheids starting 
from the measured colours. In order to test the accuracy of 
the fitted relations (Eq. [7|, we have performed a compar- 
ison of the effective temperature and gravity derived from 
the fitted relations for selected Cepheids in our sample with 
the results obtained by other authors using an independent 
spectroscopic approach. 

As an example, in Fig. |6]we have plotted the values of 
effective temperature and gravity obtained with our equa- 
tions for the pulsational cycle of the Cepheid U Sgr. In the 
same figure we also overplotted the valu es of the physical pa- 
rameters obtained spectroscopically bv lLuck fc Andrievskvl 
(2004) using the line depth ratios. We note that the curve 
of effective temperature follows the spectroscopic data with 
great accuracy. The agreement is less satisfactory for the 
curve of gravity, probably as a consequence of the cited 
problems with the W band. Similar results are obtained for 
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Figure 6. Effective temperature (bottom panel) and gravity (up- 
per panel), obtained by fitting the grid of model atmosphere, are 
plotted as function of the pulsational phase for the star U Sgr 
(P = 7.9 days). The open blue circles with the associated er- 
ror bars represent t he values of the physical pa rameters obtained 
spectroscopically bv lLuck fc Andrievskvl | |2004) . 



Figure 7. Top panel: the curve of angular diameter 6 (in mas) of 1 
Car, obtained from the surface brightness calibration using eg. 1 101 
(black poin ts), is plotted with th e interferometric m easurements 
obtained bv lKervella et al] ||2004| ) (red points) and bv lDavis et al.l 
(2009) (blue points). Bottom panel: difference between the two 
curves the mean and standard deviation of such a difference is 
labelled. 



another dozen of stars for which the comparison with spec- 
troscopic data was possible. We can conclude, that our ap- 
proximation for the effective temperature is quite accurate, 
with mean relative differences of 0.4%. On the other hand, 
we can reproduce the effective gravity, log g^ff, with mean 
relative differences of 4%. However, the larger error affecting 
the estimate of the surface gravity is not a thorny problem 
in the estimate of Cepheid radii, because the gravity enters 
only in the determination of the AB term, through the con- 
tribution of the BC to the surface brightness (Sv)- The AB 
term contributes at most with a 5% to the radius, therefore, 
we can conclude that even a significant uncertainty in our 
determination of the gravity has a minimal impact on the 
determination of Cepheid radii. 

As a further test of the accuracy of our approximation 
for the surface brightness Sv, we have compared the angu- 
lar diameter ^ of 1 Car (which is the only single Cepheid for 
which this comparison is possible) obtained through inter- 
forometric measurements (iKervella et al.l |2004| : iDavis et al.l 
[2009. ') with those obtained from the surface brightness equa- 
tion: 

e{<j}) = ioo-2(S'vW)— vW)) (10) 

where <j) is the phase along the pulsational cycle and 9 is 
given in mas (see Fig. [7]). A glance at the figure shows that 
current solution shows a systematic difference with the inter- 
ferometric measurements in the phase range between 0.2 and 
0.45. This discrepancy is due to the limited accuracy of the 
calibration we are using to estimate the surface gravity. Typ- 
ically, current calibration underestimates the gravity, in the 
above phase interval, and, in turn overestimates the radius. 
The reader interested in a thorough analysis of this problem 



is referred to IPell (|l978l ). He suggested that the main cul- 
prit of the above discrepancy is the fact that we use grid of 
atmosphere models that neglect the microturbolence vari- 
ation along the pulsation cycle. Furthermore, the problem 
in the W-band reflectivity, described above, also affects the 
surface gravity calibration. However, it is worth mentioning 
that a detailed analysis of the residuals between current es- 
timates and interferometric measurements indicate a mean 
difference and a dispersion of the order of 3%. 

Assuming that the interferometric measurements are 
minimally affected by systematic errors, we can use this av- 
erage difference to estimate our systematic indetermination 
on Sv through Eq. [T] obtaining a value of ~ 8%. Recalling 
that Sv enters only in the determination of the AB, which 
in turn is only ~ 5% of the radius, we found that the sys- 
tematic uncertainty on our estimate of Sv implies a 0.4% 
uncertainty on the radius. 

The agreement between the physical parameters ob- 
tained with our procedure and those obtained independently 
from spectroscopy and interferometry seems a valuable test 
of the accuracy in the derivation of the surface brightness 
and of the AB term, thus supporting the soundness of our 
approach. 



4 PHOTOMETRIC AND SPECTROSCOPIC 
DATA 

After the validation of the quoted approach, we applied it to 
real data in order to obtain the radii of a sample of classical 
Galactic Cepheids. 
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The photometric data and the radial velocities used in 
this work are described in the following. 

4.1 Photometric data 

The original sample we have used in this work includes 175 
Cepheids observed in the Walraven VBLUW photometric 
system during 1962 an d 1970-71 at the Leiden South ern 
Station (South Africa) l|Wahaven et al.lll964l : [Peilll976l ). It 
is 82% complete for all Cepheids brighter than V = 11.0 
magnitude at minimum of the light curve and with declina- 
tion lower than +15°. 

We have selected all the type I Cepheids (163 stars) in- 
cluded in this sample, but we have retained only those stars 
for which reliable radial velocity data were available in the 
literature (see next section). Furthermore, we have rejected 
also those stars having the colour-colour loop at least in part 
outside the grid of models. This behaviour can be observed 
for Cepheids having an extremely large loop or a loop with 
a peculiar orientation in the colour-colour plane (abnormal 
loop orientations can be caused by the presence of compan- 
ions or by peculiar chemical abundances) . A complete list of 
the selected stars, with the photometric parameters used in 
this work, is shown in Tab. [1] 

Before using the photometric data, we have corrected 
the colours f or the reddenin g effect following the procedure 
described in iPedicelli et all l|2008l ). In particular we have 
obtained the colour excess E{V — B) in the Walraven system 
from: 



E{B~V), 



2.375 -0.169(y-B) 



(11) 



E{y - B) 

where {V — B) is the mearjf] Walraven colour of the consid- 
ered star and E{B — V)j is the colour excess in the Johnson 
system. As for the colour excess in the John son filters, we 
have referred to lLanev fc Caldwelll (|2007l ') or to lFernid (|l995l ') 
when the value was not available in the former work. In 
the latter case, we have transformed the value of the colour 
excess from the Fernie's scale into the Lane y &c Caldwell 
scale b y using the relation recently provided bv lFougue et all 
(|2003), namely E{B - V')£&c = 0.9525(5 - F)Fe™e 
mag. The {U — W) colour ex cesses follow from th e relation 
E{U - W) = 0A5E{V - B) l|Pedicelli et al.ll2008l ). 



4.2 Radial velocity data 

We have obtained the radial velocity data from different 
sources available in the literature. The references of radial 
velocity data are listed in Tab.[l]for each star of the sample. 



The main catalogues used in this work are from: 


IStibbd (Il955l): GierenI (Il977 


. 19811. 19851): Coulson et all 


(Il985l): iBarnes et al. (Il987l. 


19881): iMetzEer et all (1 19931): 


iBersier et all 1 


19941): iTavlor et al. lll997f): Gorvnva et all 


I1998t): llmbert 


(19991): iBersieil (20021): IStormI (2004). The 



sample of radial velocity measurements was selected accord- 
ing to the following criteria. We neglected Cepheids with less 
than 10 measurements. For the Cepheids with multiple sam- 
ples of radial velocity measurements, we selected the data 

^ The mean is derived by applying the integral mean value the- 
orem to the colour curve. 



set characterized either by the most accurate measurements 
-typically better than ~1 km/s- and/or closer in time with 
photometric data. The latter criterion was adopted to mini- 
mize possible problems due to period changes and/or phase 
of maxima that are not very precise. In a few cases (nine 
stars), we have combined two or more data sets in order to 
improve the sampling along the radial velocity curve. To ac- 
complish this goal, we corrected for possible shifts in radial 
velocity among the different data sets. The shift was deter- 
mined by interpolating the most accurate sample of radial 
velocities and then by estimating the median distance -along 
the velocity axis- between the fitted curve and the other ra- 
dial velocity samples. The values of the shift we found are 
listed in Tab. [T] 

In order to obtain an accurate phasing between the pho- 
tometric and spectroscopic data, we recalculated the phases 
of radial velocity measurements by using the same epoch 
and period as for the photometry. The final sample used in 
this work consists of 63 Cepheids whose main characteristics 
are listed in Tab. [1] 



5 DERIVATION OF THE RADIUS OF 
CEPHEIDS 

5.1 The projection factor p 

Using the photometric data and the radial velocities dis- 
cussed in the previous section, we are able to derive the 
mean radius for each star of our sample following the pro- 
cedure introduced in Sees. [5] and [31 However, before pro- 
ceeding we have to fix the value of the projection parameter 
p. In the recent literature there is a strong debate about 
the correct value to use for this important parameter. More- 
over, it is not clear yet if p is constant over the entire period 
range covered by the Cepheids . A d etailed review of this 
subject can be found in lBarne i (|2009l ) and there is no need 
to repeat here the discussion outlined in that paper, but the 
conclusion that can be drawn is that the p value is uncer- 
tain at a level of 5-10%. The same uncertainties apply to the 
possible dependence of p on the Cepheid pulsation period. 
As a matter of fact, several au thors find either a m ild (see 
e.g. Nardetto et al. I l2007l . l2009l ) or a strong (see e.g. 'Gieren| 



200^ dependence on period. On the contrar y. ICro cncwegen 
j2007), using Cepheids with HST parallaxes (jBenedict et al. 
l2007i) and interferometrically measured angular diameters, 
find, convincingly, a value of p = 1.27 ± 0.05 with no de- 
pendence on the period. The same result (p — 1.27 ± 0.06) 
was found geometrically for S Cephei by using interferomet- 
ric measurements, obtained with the C HARA Array, and 
published parallax (|Merand et al.|[2005l ). 

Moreover, the projection factor p is usually assumed to 
be constant along the pulsation cycle. However, this assump- 
tion was questioned by Rip epi et a l. (1997), since byusing 
the changes of p predicted bv ISasselov fc Karovskai (|l994 ). 
they found that the radius determinations change by ap- 
proximately 6%. The change of the p factor along the cycle 
has b een the crossroad o f several theoretical (jMar engo et al.1 
l2002i ) and empirical ((Nardetto et al.1 12004| ) investigations, 
but we still lack firm quantitative estimates of the changes. 
This is the reason why it is neglected in recent estimates of 
the Cepheid radii. We also neglected this change to perform 
comparison with similar data available in the literature. 
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Table 1. The names of Cepheids used in this work together with their period are hsted in the first and second columns respective! 



Lanev fc Caldwell ( 2007|) 



m 



The third and fo urth co l umns contain respectively the colour excess E{B — V) in the Johnson system, from 
(L) or those from lFernid l |l995 !) (F) transformed in the Laney and Caldwell system by using the relation from 
the E{V — B) in the Walraven system calculated from Eq. I lf ID . In the last column we list the source of the radial velocity data and, 
when more than one source is taken into account the shift in radial velocity (km s^^), if any, is given in parentheses. 



Fougue et al.l (120071 ). and 



Star Period (days) E(B-V) (Johnson) E(V-B) (Walraven) Radial velocity 



7J Acjl 


1 .1 1 OD 


n 1 

U. IZD 






PM Aril 


6 1142 


u.ooo 




U.ZoO 


FN Aql 


9.4822 


0.491 


CT,"l 


0.214 


TT Aql 


13.7546 


0.417 




0.217 


U Aql 


7.0239 


0.351 




0.173 


\/'AQfk Aril 

V ^tyu .ri.qi 


U.OU 1 u 


u.ooo 




n 1 fi7 


1X1 K^rvicL 


A fiyeo 










1 9 ^^(^90 
-LZ . OUZU 


u.oyu 




U.ZOO 


AO Par 


y . i uyu 






n 07 

U.U i 


FR Par 


7 71 R7 


ft riQfi 


) 


0A9 


1 Car 


oo.ooou 


U. -LOO 


(1 "1 


U.U 1 o 


U Car 


OO . i uo ± 


0.260 


y^) 


0.124 


V Car 


6.6967 


0.170 




0.075 


VY Car 


1 S Q91 

J-O. yz j_o 






1 Ofi 
U. lUU 


yvyv KyOil 


15.7162 


0.361 




0.157 


XY Car 


12.4348 


0.416 




0.182 


A 7, Cpn 


3.2107 


0.159 




0.068 


V Cen 


5.4939 


0.305 


CT.I 


0.125 




9.4672 


0.415 




0.181 




-Lu. yofjo 


0.267 




0.116 


AG Cru 


3.8373 


0.171 


\'^) 


0.074 




^.uyuu 


144 




0fi9 

u.uuz 


X Cru 


9900 


0.272 


) 


0.118 


W Cpm 


7.9141 


0.253 




0.123 


QV Mnn 




0.220 


(1,') 


0.108 


l\, 1\1 Ho 


7.5099 


0.114 


('^ ) 


0.049 


\J KJ 1\1 U.O 


11.6364 


0.414 




0.180 


S Nor 


9.7549 


0.180 


CT,"l 


0.078 


U Nor 


12.6413 


0.816 




0.396 


RF Onh 


4.0678 


0.223 




0.096 


Y Oph 


17.1241 


0.660 


(I,') 


0.290 


RS Ori 


7.5668 


0.363 




0.169 


AP Pup 


5.0843 


0.198 


("F"! 
) 


0.086 


AT Pup 




n 907 


(1 \ 


OSQ 

u.uoy 


RQ Pun 


41.3876 


0.454 




0.197 


X Pup 


zo. yu-i-U 


404 


11 \ 


1 77 

U. J- ( ( 


WX Pun 

VV -/V 1 LIU 


8.9382 


0.303 




0.138 


J.VV OL-W 


6.0613 


0.365 


\'^) 


0.158 


V482 Sco 


4.5279 


0.334 




0.145 


V500 Sco 


9.3166 


0.618 


(L) 


0.263 


V636 Sco 


6.7966 


0.207 


(F) 


0.089 


EV Set 


3.0910 


0.703 


(L) 


0.300 


SS Set 


3.6712 


0.324 


(L) 


0.140 


X Set 


4.1981 


0.589 


(F) 


0.257 


Y Set 


10.3415 


0.751 


(L) 


0.364 


AP Sgr 


5.0579 


0.179 


(L) 


0.077 


BB Sgr 


6.6370 


0.302 


(L) 


0.131 


U Sgr 


6.7449 


0.408 


(L) 


0.176 


X Sgr 


7.0122 


0.284 


(L) 


0.122 


Y Sgr 


5.7733 


0.187 


(L) 


0.089 


YZ Sgr 


9.5534 


0.286 


(L) 


0.127 


W Sgr 


7.5947 


0.103 


(L) 


0.048 


R TrA 


3.3893 


0.144 


(L) 


0.054 


S TrA 


6.3234 


0.081 


(L) 


0.035 


AH Vel 


4.2271 


0.070 


(F) 


0.030 


AX Vel 


2.5929 


0.239 


(L) 


0.103 


BG Vel 


6.9236 


0.426 


(F) 


0.186 


RY Vel 


28.1270 


0.545 


(L) 


0.240 



Barnes et al. (1987) 
Gorvnva et al. (19981 
Gorynya et aL (1998) 
a^rsicr (20^) 
Bersier (2002) 
Gieren (1981) 

Barnes et al. Cl988) f+0.5). Gorvnva et al. (1998) 
Coulson et al. (1985) 
Coulson et al. (1985) 

Lloyd Evans (1980) (+0.4). Stibbs (1955) 
Tavlor et al. (1997) 
a^ulson et aL (1985) 
Stibbs (1955) 
Coulson et al. (1985) 
Coulson et al. (1985) 
Coulson et al. (1985) 
Gieren (1981) 
Bersier (2002) 
Coulson et al. (1985) 
Coulson et al. (1985) 
Gieren (1981)(+4.2). Stibbs (1955) 
&ercn (1981) 
_Bersier (2002) 
Imbert (1999) 
Imbert (1999) 
Llovd Evans (1980) 
Bersier (2002) 
Bersier et al. (1994) 
Coulson et al. (1985) 
_Gieren (1981) 
Gorynya et aL (1998 ) 

Gorynya et aL (1998) . Imbert (1999)(-0.3) 
Itibbs (1955) 
Gieren (1985) 
Storm (2004) 
Barnes et al. (1988) 

Stibbs (1955) . Barnes et al. (1988)(-1.7) 

Llovd Evans (1980) (+2.9). Gieren (1981) , Stibbs (1955) (-1.1) 
ffieren (1981) 
_Stibbs (1955) 
Stibbs (1955) 

Metzger e t al. (1993) (+0.4). Bersier et al. (1994) (+0.4). Storm (2004) 

Gieren (1981) 

Metzger et al. (1993) 

Barnes et al. (1988) 

Gieren (1981) 

Gieren (1981) 

Bersier et al. (1994) 

Wilson et al. (1989) 

Bersier (200^) 

Uoyd Evans (1980) : Barnes et al. (1988) 
Bersier et al. ( 1994) 
Gieren (1981) 

Llovd Evans (1980) (+1.4). Gieren (1981) 
Gieren (1977) 
Stibbs (1955) 
Stibbs (1955) 
Coulson et al. (1985) 
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Table 1 - continued 



Star 


Period (days) 


E(B-V) (Johnson) 


E(V-B) (Walraven) 


Radial velocity 


RZ Vel 


20.3969 


0.294 (L) 


0.129 


Coulson et al. Q985) 


SX Vel 


9.5499 


0.274 (L) 


0.108 


Stibbs ('19551 


SW Vel 


23.4744 


0.346 (L) 


0.153 


Coulson et al. (1985) 


T Vel 


4.6397 


0.289 (L) 


0.122 


Stibbs ^1955) 


V Vel 


4.3710 


0.144 (L) 


0.090 


Gieren (1985) 



Given this overall uncertainty, we decided to use two 
different values of tlie p factor: 1.36 and 1.27. Tliese two 
values were cliosen because tliey include tlie range of val- 
ues found in the literature and, in particular the value 1.36 
was selected to allow comparison with our previous works 
l|Ripepi et al.ll 19971 : IRuoppo et al]|2004l ). 

5.2 Radius calculation 

The FORTRAN77 code used in this work performs a fit of 
the magnitude curve, the [V — B), {U — W) colour curves and 
the radial velocity curve using a Fourier fit, with a number 
of harmonics fixed interactively by the user. Then it solves 
the CORS Eq. ([2} for the radius at an arbitrary phase both 
with and without the A_B term. The mean radius is, finally, 
calculated by integrating the radial velocity curve twice. 

Table |3] gives the results of our procedure together 
with the radii obtained by other authors. For three stars 
the CORS procedure does not converge and consequently 
the radius is not available. In particular we have listed 
the name of the star, information on binarity, the period, 
the radius obtained with and withou t the Ai? term and 
some results found in the literature (Ruoppo ct al. 2004; 
iBarnes et all l2005l : iRoio AreUano fc Arellano Fcrr o 1994). 
In order to perform the comparison with similar estimates 
available in the literature, we rescaled their radius determi- 
nations according to our value of the projection factor. We 
found that our rad ius determinations agree quite w e ll wit h 
those provided by iRoio Arellano fc Arellano Ferrol l|l994h . 
with a median difference of less than 1% and a dispersion 
of 15%. On the other hand, current r adius evaluations are 
~ 13% smaller than those provid ed bvlRuoppq et EiL ()2004l ) 
and ~ 21% larger than those bv lBarnes et al.l l|200^ ), with 
a dispersion of the order of 14% for both of them. 

As further test, we have selected the Cepheids of our 
sample with available interferometric measurements of the 
radius. In par t icular , we compared the radii contained in 
iKervella et al.l (|2004l ). obtained with p — 1.36, with our re- 
sults corresponding to the same projection factor. The stars 
in common with their sample are: X Sgr, rj Aql, W Sgr, Y 
Oph, 1 Car. As shown in Tab. [2] our radii are completely 
consistent with interferometric measurements, except for Y 
Oph, but the error given by Kervella et al. is extremely large. 

We found another estim atio n of the radius of 1 Car and 
T) Aql in iDavis et al.] l|2009l ) and IJacobI l|2008h respectively. 
Using interferometric observations made with the Sydney 
University Stellar Interferometer, they found R = 169±8i?Q 
for 1 Car and i? = 39 ± 6-Rq for rj Aql, adopting a constant 
projection factor p = 1.30. We rescaled their estimations 
using our projection factor and obtained R = 177 ± 9,Rq 
for 1 Car, which is in agreement with our result within 0.9 



Table 2. Interferometric radii derived bv IKervella et aL I ||2004| ) 
compared with the results of this work for p = 1.36. In the first 
column there is the name of the Cepheids in common with the 
sample of Kervella et al, while in the second and third columns 
there are the interferometric radius and our estimation respec- 
tively. 



Star 



Kervella et al. 



This work 



XSgr 


52.2+j^ 


52.4+1^ 


Tj Aql 




59.0l« 


W Sgr 


56.4l?« 


53.4l*-l 


Y Oph 






1 Car 


191 2+^'^ 


190. 7l}^;^ 



a, while we are not in agreement with the radius of rj Aql, 
7? = 41±6i?0. 



6 THE PERIOD-RADIUS RELATION 

Since the radii obtained in this work are consistent with 
those available in the literature, we have investigated the 
Period-Radius relation for the present data set. 

In order to derive the Period-Radius relation, we ne- 
glected the first overtone pulsators, since these objects have, 
at fixed period, systematically larger radii than fundamental 
pulsators. To pin point first overtone pu lsators i n our sample 
we followed the classification given bv lFernid (|l995l ). Data 
plotted in Fig. 8 show the radii for the entire Cepheid sam- 
ple. The different colors mark Cepheids in different types 
of binary systems. Following the classification suggested by 
ISzabadod (|2003l ). the black c rcles display single stars and 
visual binaries (V), the red points show the spectroscopic 
binaries (B), the blue points the spectroscopic binaries with 
known orbits (O) and the magenta points mark the bina- 
ries with photometric companions, but their binarity has to 
be confirmed (b), while the green circles show the position 
of the first overtone pulsators. We ended up with a sam- 
ple of 56 Cepheids and they have been used to derive a new 
Period-Radius relation for Galactic Cepheids. We performed 
a linear fit in the form log R = a log P -\-h, with the radius 
in solar units. Finally, by using p — 1.36, we have found the 
following relations: 



logR = (0.69 ± 0.03) log P + (1.17 ± 0.03) 
logi? = (0.73 ± 0.03) log P + (1.14 ± 0.03) 



(12) 
(13) 



in the case without and with the AB term, respectively. The 
standard deviation of the above Period-Radius relations is 
0.06 dex. A glance at the data plotted in Fig. [8] indicates 
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Table 3. R adii derived in the present work compared with those previously determined by other authors ||Ruopdo et al .l2004l : lBarne7 et alj 
2OO5I: iRoio Arellano & Arellano Fcrro 1994). The third column gives information on binarity according to the nomenclature from 
Szabadoalj2003 '): B=spectroscopic binary, b=photometric companion with no sure physical relation, 0=spectroscopic binary with known 
orbit, V= visual binary. The Ccpheids with no binarity label are bona fide single stars. 
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61.0 
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AP Sgr 
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48.6 


47.9 
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44.7 


62.3 




46.1 


BB Sgr 


6.6370 


B 


44.2 


43.8 


41.3 


40.9 


61.3 


45.6 


54.1 


U Sgr 


6.7449 


B 


64.6 


65.7 


60.4 


61.3 


82.5 


49.1 


55.6 


W Sgr 


7.5947 





53.2 


53.4 


49.7 


49.9 


62.6 




50.8 


XSgr 


7.0122 





54.3 


52.4 


50.7 


48.9 


59.8 




44.1 
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5.7733 


B 


55.7 


55.2 


52.0 


51.5 


56.2 




48.3 


YZ Sgr 


9.5534 


B 


66.8 


69.3 


62.4 


64.7 


95.8 




69.4 


R TrA 


3.3893 


B 


35.0 


35.1 


32.7 


32.8 








S TrA 


6.3234 




54.1 


54.8 


50.5 


51.2 








AH Vel 


4.2271 


B 


55.9 


51.7 


52.2 


48.3 








AX Vel 


2.5928 




53.8 


54.1 


50.2 


50.5 








BG Vel 


6.9236 


b 


66.0 


69.7 


61.6 


65.1 








RY Vel 


28.1270 




172.3 


181.2 


160.9 


169.2 




136.4 
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Table 3 - continued 



Star 


Period(days) 


Binarity 


r^p = 1.36 
withoutAB 


7?P = l-36 


„p = 1.27 

withoutAB 


R-AB 


R-RuoppoOA RBarnesQ5 RRojodA 






«0 




Rq Rq Rq 


RZ Vel 


20.3969 




130.8 


137.0 


122.1 


128.0 


110.6 


SX Vel 


9.5499 




64.3 


65.1 


60.1 


60.8 




SW Vel 


23.4744 




140.1 


151.9 


130.3 


141.8 


122.2 


T Vel 


4.6397 


B 


57.2 


57.3 


53.4 


53.5 


40.0 


V Vel 


4.3710 


B 


34.4 


34.2 


32.1 


31.9 





that a significant fraction of the scatter is due to binaries. 
The flux of the companions affect the radius estimates. This 
effect is more significant in optical and in ultraviolet bands 
than in near infrared bands, because t he Cepheid com pan- 
ions are typically blue stars (see e.g. ISzabadoi 120031 '). In 
order to provide a tighter Period-Radius relation, we ne- 
glected all the spectroscopic binaries from the fit (types B 
and O of the above classification) and the first overtone can- 
didate X Set (see Sec. 16. If) . The final sample reduces to 26 
Cepheids and the corresponding Period-Radius relations for 
p = 1.36 are the following: 

logi?= (0.71 ±0.03) log P 4- (1.16 ±0.03) (14) 

log R = (0.75 ± 0.03) log P + (1.13 ± 0.03) (15) 

with an rms of 0.04 and 0.03 dex in the case without and 
with the AB term, respectively. The exclusion of the bi- 
nary stars does not affect significantly the coefficients of 
the Period-Radius relation, but it significantly decreases the 
scatter of the linear fit. 

The fitted equations for the case p — 1.27 are the fol- 
lowing: 

log R = (0.71 ± 0.03) log P ± (1.13 ± 0.03) (16) 

log R = (0.75 ± 0.03) log P + (1.10 ± 0.03) (17) 

for the case without and with the AB term respectively. The 
rms errors of these relation are the same as for the case with 
p=1.36. 

In order to validate current Period-Radius relations, 
we performed a detailed comparison with similar theoret- 
ical and empirical relations available in the literature. The 
zero-points and the slopes are listed in Tab. [l] 

In order to compare theory and observations we se- 
lected three different predictions concerning the Period- 
Radius relations at solar chemical compositi on. In partic- 
ular t he Period-Radius relation provided by lAlibert et al.l 
l|l999l l by using linear, nonadiabatic, convective models con- 
structed assuming a canonical (i.e. based on evolutionary 
models that neglect convective core overshooting during cen- 
tral hydrogen burning phase, rotation and mass loss) Mass- 
Luminosity relation. Moreover , we also adopted pred i ction s 
provided bv lBono et al] (|l998l ) and bv lPetroni et all (|2003l ) 
by using, nonlinear, convective models constructed using 
both canonical and noncanonical Mass-Luminosity relations 
(mimicking the increase in luminosity caused by overshoot- 
ing). To provide a robust comparison between theory and 
observations, we performed a new estimate of the Period- 
Radius relations using only the models of the last two data 
sets covering the same period range as of the current sample 
of Galactic Cepheids. The zero-points and the slopes listed 
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Figure 8. Period— Radius relations obtained in the case with AB 
term (bottom panel) and without AB term (upper panel), using 
all the Cepheids of our sample, with the exception of first over- 
tone pulsators (green points). Single stars and visual binaries are 
plotted as black circles, spectroscopic binaries as red points, spec- 
troscopic binaries with known orbit as blue points and binaries 
with photometric companion with unknown physical relations as 
magenta points. 



in the first three entries of Tab. [3] indicate that they agree 
quite well with empirical Period-Radius relations obtained 
by neglecting the AB term. Theoretical and empirical un- 
certainty do not allow us to constrain the precision of the 
Period-Radius relations based on the two different projec- 
tion factors. 

The values given in Tab.|4]indicate that current Period- 
Radius relations agree, within the errors, with similar em- 
pirical estimates. On the other h a nd, t he Period-Radius 
relation provided by iRipepi et all (|l997l ), using the same 
method, gives a smaller zero-point, but a steeper slope. 

It is interesting to note that the zero-points of our 
relations seem to change in opposite directions by includ- 
ing or exclud i ng th e AB term with respect to those by 
IRuoppo et al] (|2004| ). The main result derived from Tab. |4] 
is that the agreement of the zero-point of our relations with 
those of other authors (excluding the theoretical result) im- 
proves by including the AB term and by decreasing the pro- 
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Table 4. Coefficients of the Period-Radius relation from tfie literature and of our work. The first two columns contain the slope and the 
intercept of the linear relation, the third column contains the reference, the fourth column the technique used to derive the Period-Radius 
relation and the last one the number of stars used in the fit. 



a 


b 


Source 


Method 


N 


0.71 


1.14 


Alibert et al. Qgog")" 


Theory (Canonical models; Solar metallicity) 




0.673 ± 0.009 


1.164 ±0.010 


Bono et al. C 1998) 


Theory (Canonical+Noncanonical models; Solar metallicity) 




0.699 ±0.011 


1.150 ±0.012 


Petroni et al. f2003") 


Theory (Canonical+Noncanonical models; Solar metallicity) 




0.75 ±0.02 


1.07 ±0.02 


Gieren et al. fl998')'' 


Surf. Brightness (variable p) 


28 


0.74 ±0.03 


1.12 ± 0.03 


Roio Arellano & Arellano Ferro (1994) 


Surf. Brightness (p = 1.31) 


54 


0.606 ± 0.037 


1.263 ±0.033 


Ripepi et ah (1997) 


Surf. Brightness (AS ^ 0, p = 1.36) 


64 


0.767 ±0.009 


1.091 ±0.011 


Kervella et al. (2004) 


Interferometry (p = 1.36) 


8 


0.686 ± 0.036 


1.134 ±0.034 


Grocnewegen f2007) 


Interferometry (p = 1.27) 


5 


0.747 ±0.028 


1.071 ±0.025 


Turner & Burke (2002) 


Modified Baade-Wesselink (p = 1.31) 


13 


0.69 ± 0.09 


1.18 ± 0.08 


Ruoppo et al. (2004) 


New CORS (AB = 0, p = 1.36) 


20 


0.74 ±0.11 


1.19 ± 0.09 


RuoDDO et al. (2004) 


New CORS (AB 7^ 0, p = 1.36) 


16 


0.71 ±0.03 


1.16 ± 0.03 


This Work 


New CORS (AB = 0, p = 1.36 ) 


26 


0.75 ±0.03 


1.13 ± 0.03 


This Work 


New CORS (AB 7^ 0, p = 1.36) 




0.71 ±0.03 


1.13 ±0.03 


This Work 


New CORS (AS = 0, p = 1.27 ) 




0.75 ±0.03 


1.10 ±0.03 


This Work 


New CORS (AS 7^ 0, p = 1.27) 





° Errors on the coefficients are not indicated in the original paper j Alibert et al ] l|l999l ). Tab. 5). 

'in order to have uniformity in the comparison with other relations, the listed uncertainties are the standard deviations on the 
parameters, which have been recalculated by fitting their data. 




Figure 9. Period-Radius relations for the final sample of 
Cepheids in the case with AB (top panel) and without AB (bot- 
tom panel). In each panel the fitting equation and the scatter 
around the fit are also shown. The possible first overtone X Set 
is also plotted (filled circle) but it is excluded from the fit. 

jection factor from 1.36 to 1.27. Therefore, we will consider 
the equation: 

logi?= (0.75 ±0.03) log P+ (1.10 ±0.03) (18) 

as our best fit Period-Radius relation. A plot of our relations 
for p = 1.27 is shown in Fig. |9l the plot for p — 1.36 is not 
shown because the points are only shifted vertically by the 



Figure 10. Our best fit Period— Radius rel ation (red sol i d line ) 
is plotted together with tho se obtained b y [Ali bert et all lll999l ) 
(theor ic relation, do t line) , iGieren et al. 1 lll998l) (short dashed 
line). [R uoppo et al.l | |2004| ) (short das hed long da shed line), 
iRipepj et al.l l|l997l ) (lon g dashed line) iKervella et al.i (2004) (dot 
short dashed line) and iGroenewegenI 1I2OO7I ) (dot long dashed 
line). The error bar represents the standard deviation of the cur- 
rent relation (0.03 dex). 



same offset of 0.03. The projection factor does not affect the 
slope of the Period-Radius relation only when it is assumed 
to be constant over the period range. A visual comparison 
of Eq. [18] with other relations is shown in Fig. 1101 
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6.1 X Set: a possible first overtone pulsator 

To test the first overtone nature of X Set, we estimated 
tlie radius predicted from our Period-Radius relations for a 
Cepheid with the same period: R — SG.OljX for P ~ 1-27 
and R = SO.Glj'g with p = 1.36, in units solar radii. The 
values obtained from the CORS method are R = 50.0 and 
R = 53.5, forp = 1.27 andp — 1.36 respectively (see Tab. [3]). 
They differ by more than 4(t from the predicted values, being 
larger, as expected for a first overtone pulsator. 

As a further test we have also estimated the radius 
predicted by the theoretical Pe riod-Radius relation for first 
overtone Cepheids obtained bv lBono et al.l (|2001al ): 

log R = (0.755 ± 0.007) log P + (1.250 ± 0.005) (19) 
log 7? = (0.737 ± 0.005) log P + (1.219 ± 0.004) (20) 

for canonical and non-canonical models respectively. The ra- 
dius corresponding to the period of X Set is _R = 52.5't\''[ 
and R = 47. 7^^;^ using the canonical and the non-canonical 
relation respectively. It is evident that the value obtained in 
the case of the canonical relation is in excellent agreement 
with that obtained from CORS with p = 1.36. The value ob- 
tained from non-canonical relation is consistent only within 
- 3.3(T with the CORS value for p = 1.27. 

Furthermore, we have also compared the radii of X Set 
with those of the first overtone Cepheid AH Vel (P = 4.23) 
and of the fundamental pulsator V Vel (P = 4.37), obtained 
using CORS (see Tab. [3]). As a result we find that the ra- 
dius of X Set is much more consistent with a first overtone 
pulsation. 

On the basis of this test, we are confident that X Set 
can be a good first overtone candidate, but further analysis 
is required to confirm this hypothesis. 



7 THE PERIOD-LUMINOSITY RELATION 

The estimated radii and effective temperatures are used to 
derive the intrinsic luminosities of the investigated Cepheids 
and, in turn their absolute visual magnitude. To this aim 
we assum ed the solar visual magnitude in the V band to 
be -26.75 (|Havesll 19851 ). which gives an absolute visual mag- 
nitude of -1-4.82 and an absolute bolometric magnitude of 
+4.75, with a bolometric correction of -0.07. Furthermore, 
we calculated the bolometric correction from Eq.|Xl In order 
to evaluate the absolute magnitude uncertainty, we assumed 
the scatter around the fitted logre//(y — B,U — W) and 
BC{V — B,U — W) surfaces as the uncertainty of the ef- 
fective temperature and bolometric correction respectively. 
The values of absolute magnitudes and their uncertainties, 
SMv, are listed in Tab.[5]for the case p — 1.27. 

The Period-Luminosity relation is obtained by a 
weighted linear fit of the relation Alv ~ a log P + b and 
the resulting coefficients are listed in Tab. [6] together with 
some results from other authors. We list the coefficients ob- 
tained using both values of the projection factor, p — 1.27 
and p = 1.36. 

As for the Period-Radius relation, the projection fac- 
tor influences only the zero-point of the Period-Luminosity 
relation, leaving the slope unchanged. 

From Tab. [6] it is evident that the slope of our Period- 
Luminosity relation is steeper than those of the relations 
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Figure 11. Bottom panel: our Period-Luminosity relation 
(p = 1.2 7) (red solid l i ne) is plotted togeth er with those ob - 
tained bv lKervella et al.1 1|2004|) fdo tted line), B arnes et all l|2003l ) 
(short dashed line) , IFouquT^ al.l ll2007l ) (long dashed line) and 
iGieren et al.l l|l998h (dot short dashed line). Upper panel: the 
Period-Wesenheit relation obtained using the projection fac- 
tor p = 1.27 (red solid line). The othe r two plotted fits are 
from lFououe et al] ||2007| ) (dotted line) and lBenedict et aL 1 1I2OO7I ) 
(short dashed line). The rms around the fitted relations are indi- 
cated in both panels. 



bv iFouaue et al.l (|2007l ) and iBarnes et al.l (|2003l ). although 
they are consistent within the errors. On the other hand, it 
is in excellen t agre emen t with the slop e of th e relations by 
iGieren et aO (| 19981 ) and lKervella et al.1 (j2004j ). 

The agreement of the zero-point with other works is bet- 
ter if we choose the projection factor p — 1.27, while in the 
case with p — 1.36 the Period-Luminosity relation seems 
to make Cepheids too luminous, par ticularly if w e com - 
pa re our relation with that found by iFougue et al.l (|2007l ) 
or lCieren et al.l (|l998l ). We note that the systematic errors 
on the zero-point of the Period-Luminosity relations ob- 
tained for the two values of p is ~ 0.15 mag, i.e. comparable 
with the statistical error on the zero-point itself. 

A plot showing the comparison of our Period- 
Lumin o sity relation (p = 1.27) w ith those fromfearnes et al.l 
ll2003l). iKervella et al.l (j2004h . iFouaue et all l|2007l ) and 
ICieren et al.l l| 19981 ) is represented in Fig. 1111 

Using the V band absolute magnitudes and 
intensity- mea n ((V ) — (Ic)) colour tabulated 
Grocnewe genI l|l999t ) for all our Cepheids, except 
V500 Sco, W Gem, WX Pup, X Cru and Y Set, 
which we used the values from ICaldwell fc CoulsonI 
(|l987h , we derived the reddening free Wesenheit mag- 
nitude Wvic = Mv - RviciiV) - (Ic)), where 
the colour coefficient Rvi r- = 2.55 is adopted from 
ICaputo. Marconi fc Musellal (|200G| ) (see Tab. [5] for Wvic 
values and their uncertainty SWvic)- The Period-Wesenheit 
relation, Wvic ~ alogP + b, has been derived through 
a weighted linear fit using both values of the projection 
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Table 5. The data used to fit the Period— Luminosity and Period- Wesenheit relations are hsted for the Cepheids of our sample. The 
absolute magnitudes My u sed to derive t h e Wes enheit magnitudes are calculated as described in the text, while the {{V) — (Ic)) colour 
in the last column i s from 'G roenewegenl lll999l 'l. For V500 Sco, W Gem, WX Pup, X Cru and Y Set the {{V) — (Ic)) colour is not 
available in lGroeiiewegem il999l ') and we have used the value from iCaldwell fc CoulsonI l ll987t) . 



Star name 


Period (days) 


'^eff (K) 


My (mag) 


5Mv (mag) 


Wvic (mag) 




{{y) - Uc)) (mag) 


FM Aql 


6.11423 


5746 


-3.50 


0.17 


-5.60 


0.18 


1.531 


FN Aql 


9.48224 


5731 


-4.22 


0.18 


-6.24 


0.19 


1.389 


TT Aql 


13.7546 


5293 


-4.35 


0.19 


-6.65 


0.19 


1.409 


AQ Car 


9.76896 


5627 


-4.15 


0.18 


-6.12 


0.19 


0.982 


ER Car 


7.7187 


5642 


-3.78 


0.20 


-5.69 


0.21 


0.865 


1 Car 


35.5330 


4843 


-5.35 


0.22 


-7.93 


0.22 


1.177 


XY Car 


12.43483 


5503 


-4.41 


0.19 


-6.55 


0.19 


1.342 


V Cen 


5.49392 


6069 


-3.62 


0.17 


-5.33 


0.17 


1.041 


V339 Cen 


9.4672 


5591 


-4.07 


0.18 


-6.13 


0.19 


1.314 


S Cru 


4.68997 


5956 


-3.27 


0.17 


-4.98 


0.17 


0.846 


X Cru 


6.219970 


5716 


-3.50 


0.19 


-5.44 


0.20 


1.090 


W Gem 


7.914130 


5854 


-4.03 


0.18 


-5.90 


0.18 


1.043 


SV Mon 


15.2321 


5391 


-4.63 


0.19 


-6.80 


0.20 


1.119 


UU Mus 


11.63641 


5648 


-4.46 


0.19 


-6.47 


0.19 


1.294 


U Nor 


12.64133 


5476 


-4.41 


0.19 


-6.66 


0.19 


1.874 


RS Pup 


41.3876 


5087 


-5.91 


0.22 


-8.41 


0.23 


1.531 


X Pup 


25.9610 


5368 


-5.46 


0.21 


-7.69 


0.21 


1.365 


WX Pup 


8.938250 


5809 


-4.18 


0.18 


-6.07 


0.18 


1.108 


V500 Sco 


9.316650 


5828 


-4.28 


0.18 


-6.25 


0.19 


1.527 


SS Set 


3.671253 


5970 


-2.88 


0.16 


-4.65 


0.17 


1.086 


Y Set 


10.341504 


5523 


-4.14 


0.18 


-6.43 


0.19 


1.810 


S TrA 


6.32344 


5810 


-3.62 


0.17 


-5.40 


0.18 


0.796 


RY Vel 


28.1270 


5411 


-5.65 


0.21 


-7.86 


0.21 


1.532 


RZ Vel 


20.3969 


5316 


-5.01 


0.20 


-7.18 


0.21 


1.208 


SX Vel 


9.54993 


5992 


-4.46 


0.21 


-6.15 


0.22 


0.996 


SW Vel 


23.4744 


5355 


-5.29 


0.21 


-7.46 


0.21 


1.272 



Table 6. Coefficients of the Period-Luminosity and of the Period- Wesenheit relations obtained in this work and from other authors. 



Band 


a 


b 


rms 


Source 


V 


-2.77 ±0.07 


-1.29 ±0.08 


0.20 


(Gieren et al. 1998) 


V 


-2.77 ±0.07 


-1.44 ±0.10 




(Kervella et al. 2004) 


V 


-2.68 ±0.08 


-1.27 ±0.02 


0.17 


(FouQue et al. 2007) 


V 


-2.69 ±0.17 


-1.47 ±0.28 


0.19 


(Barnes et al. 2003) 


V 


-2.78 ±0.11 


-1.57 ±0.11 


0.13 


This Work (p=1.36) 


V 


-2.78 ±0.11 


-1.42 ±0.11 


0.13 


This Work (p=1.27) 


Wvic 


-3.48 ± 0.07 


-2.41 ±0.02 


0.17 


(Fouque et al. 2007) 


Wvic 


-3.34 ±0.17 


-2.52 ±0.17 


0.09 


(Benedict et al. 2007) 


Wvic 


-3.50 ±0.07 


-2.87 ±0.07 


0.09 


This Work (p=1.36) 


Wvia 


-3.50 ±0.07 


-2.72 ±0.07 


0.09 


This Work (p=1.27) 



factor; the resul ts are listed i n Tab. [6] toge ther with those 
obtain ed from iFougue et al.l (120071 1 and [Benedict et al.l 
l|2007tl . 

The slope of our Wesenheit relation is in agreement 
with the results from other authors, using both values of 
the projection factor. The variation of the projection fac- 
tor influences the zero-point, which, for p — 1.36 seems to 
be too luminous , up t o 0.46 mag if we consider the one 
by iFouque et al.l (|2007l ). This result, obviously, reflects the 
behaviour of the obtained Period-Luminosity relation (see 
above). For p = 1.27, the Wesenheit z ero-point results to be 
brighter by 0.20 mag than the one by [Benedict et all (|2007| ) 
and by 0.31 mag than the one bv lFougui^ et all (|2007l ). 



The implications of these results for the extragalactic 
distance scale calibration will be addressed in a forthcoming 
paper (Molinaro et al. in preparation). 



8 CONCLUSIONS 

In this work we have developed a modified version of the 
CORS method in the Walr aven filters, fol l owing the ap- 
proach already adopted by IRuoppo" et al.l l|2004h in the 
Stromgren system, to derive the radii of a sample of 63 
Galactic Cepheids observed in the Southern hemisphere. 
The adopted photometric data are the results of an exten- 
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sive observationa l program (|Walraven et al.lll964l : [Pelll 19761: 
iPel fc LubllioOTi '). whereas for the radial velocities we have 
used data available in the literature from various authors. 

The observations were performed in the Walraven five 
bands system (V, B, L, U, W) well-suited to measure spec- 
tral characteristics of stars like Cepheids and RR Lyrae. 
Even though the Walraven photometric system is not avail- 
able in modern observational facilities, the possibility to use, 
in the near future, the tunable filters opens a new path to- 
ward the performance of Cepheid observations in the Wal- 
raven bands. Since new stellar atmosphere model colours 
in the Walrav en photometric system were available from 
ICasteUil (1999), we have been able to calibrate the surface 
brightness function using {V — B) and {U ~ W) colours, in 
order to apply the complete CORS method, including the 
AB term {Eq.^. 

The comparison of temperature curves for some 
Cepheids of our sample with the temperature values ob- 
tained from spectroscopy has revealed that the (V — B) Wal- 
raven colour is a very good temperature indicator, whereas 
a worse agreement is found for the gravity curves, probably 
due to a lower accuracy in the measurement of the W band 
flux. 

In this work we are able to accurately calibrate the sur- 
face brightness and to derive the radius even for Cepheids 
with P ^ 40 days. The presence of long period stars in 
our sample allowed us to better constrain the coefficients of 
the Period-Radius relation. In order to derive the Period- 
Radius relation, we have considered only single stars, ex- 
cluding binaries, since orbital motion and/or the influence 
of the flux of a companion could afi'ect the radius derivation. 
Furthermore we have also excluded first overtone pulsators. 
Our final sample thus reduces to 27 Cepheids. 

In order to take account of the uncertainty in the pro- 
jection factor value, we have considered two possible values, 
p = 1.27 and p = 1.36, chosen in order to include the range 
of values found in the literature. 

The Period-Radius relations obtained in the present 
work are the following: 



log R = (0.75 ± 0.03) log P -f (1.10 ± 0.03) 
logi? = (0.75 ± 0.03) logP -f- (1.13 ± 0.03) 



(21) 
(22) 



for p = 1.27 and p — 1.36 respectively. We excluded X Set 
from the fit because of its large deviation (> 4a). We in- 
vestigated the hypothesis that it can be an unknown first 
overtone Cepheid and our results suggest that it is a good 
candidate, but further analysis is required. 

The comparison with independent results in the lit- 
erature has shown that the agreement with the slope of 
our Period-Radius relations improves including the AB 
term. The only exc eption is the Period-Radius relation by 
iRipepi et af] (|l997l ). who find a flatter slope using the sur- 
face brightness technique. Furthermore, the uncertainties in 
our relations are of the same order of magnitude or even 
smaller than those obtained using the infrared data. The 
zero-point of our relation decreases from b — 1.13 to 6 = 1.10 
by varying the projection factor from p = 1.36 to p — 1.27, 
providing a better agreement with several authors. There- 
fore, we consider log R = (0.75 ± 0.03) log P -|- (1.10 ± 0.03), 
obtained for p = 1.27, as our best fit Period-Radius relation. 

This result is also confirmed by the study of the 
Period-Luminosity relation. The comparison with the re- 



sults of other authors (ICieren et al. 19981: Kervella et al.l 
|2004 iFouaue et all l2007l : iBarnes et all |2003| ) shows that 
the agreement improves by decreasing the projection fac- 
tor values from 1.36 to 1.27. Thus we decided to keep the 
value p = 1.27 and considered the corresponding Period- 
Luminosity relation as the best fit one: 



Mv 



(-2.78 ± 0.11) log P - (1.42 ± 0.11) 



(23) 



with a scatter of 0.13 mag which is smaller than other works. 

In order to exclude possible uncertainty in the redden- 
ing, we also derived the Period- Wesenheit relation. As ex- 
pected from our findings about the Period-Luminosity re- 
lation, the result obtained with p = 1.36 seems to make 
Cepheids t oo luminous with r e spect to other results in the 
literature (|Fouaue et al.ll2007l : [Benedict et"al]|2007h . Thus 
we decided to keep the relation for p = 1.27 again as the 
best fit one: 



Wvic = (-3.50 ±0.07) log P- (2.72 ±0.07). 



(24) 



which has a scatter of 0.09 mag. 

Whereas the implications of the obtained relations for 
distance scale applications and for the very important issue 
of the universality of the Period-Luminosity relations de- 
serve further investigation and will be the subject of a forth- 
coming paper (Molinaro et al. in preparation), it is worth 
noticing that the presented comparison with independent 
results has provided a test for the projection factor, suggest- 
ing a smaller value than the classical p — 1.36. In particular, 
p = 1.27 appears to be the most accur ate, in agreeme nt with 
the va l ues obtained by ICroenewegenl (jl99S ) ; .Merand et al.l 
(|2005l ): lNardetto et all (|2007D . 
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APPENDIX A: THE FIT OF THEORETICAL 

GRIDS 

Here we show the result of the polynomial fit performed in 
order to express the effective temperature, log Te/ / and the 

effective gravity, log g^ff as function of the {V — B) and 
{U — W) colours and the Bolometric Corrections, BC, as 
function of temperature and gravity: 

\ogT^ff = ao+ai{V-B) + a2{U-W)+ (Al) 
+a3{V -B){U -W) 

logffe// = ho + hi{V -B)+ h2(U - Wf + (A2) 

hsiV - B){U -W) + b4(V - Bf(U -W) + 
+h{V - B){U - wf + -t-66(V - Bf{U - wf 

BC = co + ci(logTe//) + C2{logT,fff + (A3) 

C3(logfl'e//) + C4(l0gTe//)(l0gSe//) 



16 R. Molinaro et al. 

Table Al. Coefficients of the polynomial relations described in the Appendix |A] 



ao 

3.959 ±0.004 


ai 

-0.666 ±0.012 


12 

-0.251 ±0.013 


"3 

0.82 ±0.03 


a4 


0-5 




bo 

3.21 ±0.19 


bi 

23.0 ±3.0 


b2 

-24.0 ±2.0 


b3 

-130.0 ± 18.0 


b4 

59.0 ± 12 


bs 

204.0 ± 27 


be 

-151.0 ±29 


CO 

85.7 ± 1.4 


Cl 

-44.41 ±0.73 


C2 

5.75 ±0.10 


-0.12 ±0.02 


C4 

0.035 ±0.006 


C5 


C6 




V-B 



Figure Bl. Grids of models in the (V — B) — (L — U) plane for 
the parameter range 5000 ^5 T'eff < 7500 and 0.0 ^ loggejj ^ 
4.00. The coloured lines represent the loop of the two long period 
Cepheids 1 Car (P = 35 days) and x Pup (P = 26 days). 



this finding in a subsequent work. A possible reason may be 
that the L band covers an extremely line-rich region of the 
spectrum. Apart from a very large number of metallic lines, 
it includes the higher Balmer lines and molecular lines of 
CN and CH. For this reason, L may be more sensitive to 
dynamic effects (emission components due to Shockwaves, 
enhanced and time-variable microturbulence) not included 
in the physics at the base of models, hence more sensitive to 
breakdown of the quasi-hydrostatic equilibrium assumption. 
The presence of molecular lines in L may also mean that 
the Kurucz/Castelli model atmospheres have shortcomings 
at the cool end of the grids that become particularly pro- 
nounced in L. 

Using only the Cepheids with loops fully contained on 
the grids, we have estimated their radii comparing them with 
old values obtained using the {U — W) colour. The ratios of 
their values, plotted in Fig. IB2I show that the differences 
are smaller than 10%, with mean ratio equal to 0.98 and 
a standard deviation of 0.03. We have also overplotted the 
new radii on the Period-Radius relation (Eq. llSp and they 
seem to lie on the fitted relation as shown in Fig. IB3I 

On the basis of these analysis, we have decided to use 
the {U — W) colour instead of {L — U). 



The coefficients ai, hi, a of the previous equations are listed 
in Tab. lAll The rms of the previous relations are 0.0033 dex, 
0.15 dex and 0.0032 respectively. 



APPENDIX B: THE ANALYSIS WITH 

{y - B)-{L - U) COLOURS. 

Figure IB II shows the grid of models in the plane {V — B) — 
[L — U) for the parameter range 5000 ^ T^ff < 7500 and 
0.0 ^ log Qeff ^ 4.00. According to the procedure described 
in Appendix [X] we performed a fit of the grid using poly- 
nomial equations in the variable (V — B) and {L — U). For 
brevity, we do not give the explicit expressions of the fit, but 
we recall that the rms of residuals around the fit are 0.0023 
and 0.14 for logTe// and log g^ff respectively. 

In the same graphic the loops of the Cepheids 1 Car 
{P = 35 days) and X Pup (P — 26 days) are also plotted. 
Since part of the loops falls outside of the grid, we are not 
able to accurately estimate the structural parameters of the 
two stars for all phases of the pulsational cycle. The same 
behaviour is observed for most Cepheids with period longer 
than 11 days. It would be interesting to investigate about 
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Figure B2. The ratios between the new radii obtained with {L — 
U) colour and the old values for a sample of singular Cepheids 
with colour-colour loop fully contained on the {V — B) — {L — 
U) grid of models. The mean ratio and the standard deviation, 
indicated in red, are respectively m = 0.98 and a = 0.03. 
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Figure B3. The Period-Radius relation (Eq. I15I I is plotted to- 
gether with the values of radius (blue stars) obtained with (L — U) 
colour, for the sample of Cepheids having loop fully contained on 
the grid of (V - B) - {L - U) models. 



